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Abstract 

We propose a new texture for the light neutrino mass matrix. The proposal is based upon imposing zero- 
trace condition on the two by two sub-matrices of the complex symmetric Majorana mass matrix in the flavor 
basis where the charged lepton mass matrix is diagonal. Restricting the mass matrix to have two traceless 
sub-matrices may be found sufficient to describe the current data. Eight out of fifteen independent possible 
cases are found to be compatible with current data. Numerical and some approximate analytical results are 
presented. 

PACS numbers: 14.60.Pq; ff.30.Hv; f4.60.St 

1 Introduction 

The observed phenomenon of neutrino oscillations [l]-[5], solar and atmospheric, provides a compelling evidence 
that neutrinos are massive and lepton flavors are mixed. These facts are in contrast to the Standard Model 
of Particle Physics, especially the Electro-weak interaction which is based on SU{2)l x U{1)y gauge, where 
the neutrinos are massless. We will assume that the neutrinos are of Majorana type neutrinos, as favored by 
some theoretical considerations [6], whence the mass matrix M is symmetric. In the frame work of three lepton 
families, the mass spectrum and flavor mixing are fully described by twelve real parameters: three charged lepton 
masses (me, m^, m^), three neutrino masses (mi, TO2, "t-s), three flavor mixing angles [OxT&yT&z), one Dirac-type 
CP-violating phase (S) and two Majorana-type CP-violating phases (p and a). 

In the flavor basis where the charged lepton mass matrix is diagonal, the mass term for Majorana neutrinos 
in terms of gauge eigen states, in the case of three flavors, has the form 

Cmass = -^'yIC-'M^i^L + h.c, (1) 

where C is the charge conjugation matrix. The complex symmetric Majorana neutrino mass matrix M can be 
diagonalized by unitary transformation that links the gauge and mass eigen states as: 

(2) 
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The mass term can be written in terms of mass eigen states as 



-Cmass = C ^ Dv' L + h.c, (3) 

where 

V ^ D (4) 

and D = diag(TOi, m2, TO3) with nii real positive numbers. The lepton flavor mixing matrix V contains six 
real parameters, three of them are mixing angles while the rest are three CP- violating phases. Following the 
parameterization in [7], the matrix V can be expressed as a product of the Dirac-type flavor mixing matrix 
U (consisting of three mixing angles and one CP-violating phase) and a diagonal matrix P (consisting of two 
nontrivial Majorana phases): V = U P. Then we may rewrite A'l in equation (4) as 

(5) 

where two Majorana- type CP-violating phases arc included into the complex neutrino mass eigenvalues Ai, and 
the relation \Xi\ — mi holds. Without loss of generality, we take 

Ai=mie^''' , A2= 7712 6^*'^ , A3 = ms (6) 

The matrix U is the Dirac-type flavor mixing matrix [7] 

U I Cx Sy Sx Cy C Sx Sy Cx Cy C Sy Cz | (7) 

^x Sz -f" Sx Sy C Sx Cy Sz Cx Sy C Cy Cz 




where we used the convention = sin a , Cq. = cos a. 

For the parameterization followed in the present work, the mixing angles {Ox, 0y,9z) are directly related to the 
angles of solar, atmospheric and CHOOZ reactor oscillations [7]: 

Ox « 6*801' ~ 6'atm' ~ ^'chz' (^) 

and 

a2 12 2|a2 12 21 /n^ 

^"^sol " r"2 - "^1 h ^"^atm = P3 - "^1 1 ■ (9) 
Beta decay, ncutrinolcss double-beta decays and precision cosmology are sensitive to the absolute neutrino 
mass scale. The dependence can be respectively characterized through two non-oscillation parameters and mass 
sum parameter as follows: the effective electron neutrino mass term 

Ml, = ^\Uei\^ml + \U,2\^ml + \U,3\^ml (10) 
the effective mass term of neutrinoless double beta decay 



M,3f3 = 7713 

and the sum mass parameter 



ms TO3 



(11) 

mi + m-z + TO3. (12) 
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A recent global analysis of neutrino oscillation data[8], at the confidence level of 95%, gives the best estimates 
of the oscillation parameters as 

Mtm = (2.4t°:^) X 10-W, 

Am2 J = (7.92 ± 0.7) X 10-5eV^ 

sin^e^ol = 0.314t°:°i; (13) 

sin'^atm = OAAt^fi, 

sin^^chz = 0.9t2-3 X 10-2, 

where ^atm ^^'^ ^sol angles relevant for atmospheric and solar neutrino oscillation respectively. While, 

^chz angle relevant to CHOOZ reactor experiment for neutrino oscillation [5]. A useful parameter i?^ can 

be defined as: 

sol 

2 

'atm 

which has possible values ranging, at the level of confidence 95%, as 



Ami 



0.025 <R^< 0.049, (15) 

this constraint constitutes a very tight criteria for the model to be phenomologically acceptable. 

Whereas the three masses of charged leptons (me, m^, m^) have precisely been measured[9], we have, concern- 
ing the absolute neutrino mass scale, only experimental bounds for Mp, Mpp and S as follow (see [8] and refs. 
therein) 

< 1.8 eV ] 

Mpp = 0.581^16 cV \ 95% confidence level. (16) 
E < 1.4 cV J 

The lower bound for disappears in case of absence of the ncutrinoless double-beta decay. At present, there 
is no available precise information on any of the CP violating phases. 

The present available data on neutrinos, even those in the foreseeable future, can not fully determine all the 
parameters in the neutrino mass matrix. A challenging theoretical task is to find out a Majorana mass matrix 
of the light neutrino consistent with the current data as summarized in equation (13). Several attempts have 
been made to obtain phenomenologicaly acceptable patterns of the neutrino mass matrix M, such as texture zeros 
[10]-[11], Zero sum condition [12], and determinant zero requirement [13] for the mass matrix. There are also 
many proposal for the mass matrix based on some symmetry group as in [14]. 

In this work we impose the condition that the trace of two possible 2x2 submatices is zero. The traceless 
condition can be considered as a non trivial generalization of the zcro-tcxturcs, since that, a zero-element can 
be viewed as a zero-trace of a 1 x 1 sub-matrix. Taking these submatriccs in pairs we obtain 15 independent 
possibilities. Out of these, we find just 8 of them to be phenomenologicaly acceptable. However there are other 
five cases which can be considered to be marginally accepted. The numerical and some analytical approximate 
results are reported for all these thirteen cases. 

The plan of the paper is as follows: in section 2, necessary formulas for the calculations are introduced beside 
a classification of the resulting mass patterns. Sections 3, 4 and 5 are respectively devoted to the resulting 
three possible mass patterns. For each pattern, we present the textures of M with two independent vanishing 
traces and compute the expressions of the two neutrino mass ratios and the Majaorana phases and other relevant 
parameters. Numerical and some approximate analytical results are presented. Consistency of models with 
experimental results are discussed. Wc end up by conclusions and discussions in section 5. 

2 Fifteen Possible texture with two Traceless submatrices 

As Af is 3 X 3 symmetric complex matrix, it totally has six independent complex entries. If we impose the 
condition that the trace of two possible 2x2 submatrices is zero, then we have in total 6 independent submatriccs 
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with zero trace. When these independent submatrices are taken into pairs, we get 15 possibihties that can be 
written as 



Mrs + M,j 
Ma/3 + M„„, 



(17) 



where each subscript runs over e(l), /i(2) and t(3), but (rs) ^ (ij) and (a/?) 7^ (nm). 
Using equation (5) wc then obtain the following constraint relations: 



J2 {Uri Usi + Ua Uji) Xi 
1=1 
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(18) 



1=1 



The solutions of equation (18) can be written as: 

Ai 

A3 
A2 
A3 

where 



ai 
k 



as h - 02 63 
h 02 - ai 62 
Qi ^3 - Q3 bi 
bi 02 — oi 62 



Url Usl + U^l Ujl 
Ual U/31 + Unl Uml- 



(19) 



(20) 



One can observe that the left-hand sides of equation (19) contain Major ana- type CP-violating phases, while the 
right-hand sides contain the Dirac-type CP-violating phase. Therefore two Majorana phases must depend upon 
the Dirac-type CP-violating phase. This dependence results simply from the texture of tracelss submatrices that 
we have taken. 

Comparing equation (19) with equation (6), we get the two neutrino mass ratios : 



mi 
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- 02 63 
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- ai 62 
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(21) 



and the two Majorana phases: 
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(22) 



With the inputs of three flavor mixing angles and the Dirac-type CP-violating phase, we would be able to predict 
the relative magnitude of three neutrino masses, the values of two Majorana phases and R^. The absolute neutrino 
mass scale can then be predicted by matching, for example, the value of Am^^j. The other remaining parameters 
Mf3, Mpf} and E can also be predicted. This predictability allows us to examine whether the chosen texture of AI 
with two traceless submatrices is empirically acceptable or not. The input values of 6^ and dy should be consistent 
with the boimds given by equation (13), typical ones can be taken as 0^ = 34'' and dy = 42''. The Dirac-type 
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Figure 1: R,^ as a function of 5, whereas 9x = 34*^, 9y = 42" and 9z = 5^^ for the model Di. For illustrative purposes the vertical 
range is restricted to the interval from zero to one. The sharp singling out of 5 = 93'' and 5 = 268" as consistent values with the 
relation 0.025 < i?^ < 0.049 is clear. 



CP-violating phase is not constrained. The strategy fohowed to obtain a good choice for 5 is to plot the parameter 
as a function of S while maintaining 9x ~ 34°, 9y = 42° and 9z ~ 5°. The constraint in equation (15) turns 

out to be generically very selective for the appropriate choice of S, as shown in Fig. 1 for a particular model Di 

(according to the nomenclature explained later). 

The resulting mass patterns turn out to be classified according to the following three classes: 

• Degenerate case which is characterized by mi ~ m2 ^ ttit, and is denoted by D. 

• Normal hierarchy which is characterized by mi ~ 7712 < and is denoted by N. 

• Inverted hierarchy which is characterized by mi ^ 7712 > and is denoted by I. 

In all our subsequent discussion we follow this order and nomenclature. 

To work out the explicit expressions of A1/A3 and A2/A3 in each case, we adopt the parameterization given 
in equation (7) for the Dirac-type flavor mixing matrix, from which we then obtain the analytical results for 
TOi/wa, 777.2 /m3, P, f, i?^, Mp and Mp/^. However the analytical results as well as the approximate one (for 
small Sz) are too lengthy to be displayed here, and thus we quote only the numerical results based on the 
exact formulae in equation (19). Wc present, for each model, the expressions of a's and 6's coefficients and the 
analytic approximate results for R^. However for the sake of presentation, the approximate analytical results of 
all parameters are presented for just one case for each pattern class. 

3 Degenerate models 

Pattern Di: Alee + ^^/i^i = 0, A/g^ + Al^^r = 0. In this pattern the required quantities a's and b's as given by 
equation (20) are 

2 2 I / -iS\2 2 2,/ I —iS\2 2 , 22 

0-1 ^ c^c^ + {-c^SySz - s^CyC ), s^c^ + (-SxSySz + aj:Cye ), a3 = s^ + SyC^, 

bl - CxC2;{~CxSyS 2; ~ SxCyG ) ~\- {~CxSySz ~ SxCyC ) ( ~ C-j;; Cy S2: SxSyC ), 

62 " ^x^zi^ '^x^y^z ~l~ C-xCyG ) H~ {j^^x^y^z C-xCyC )( SxCyS^ CxSyC ), 

bs = SzSyCz + Syclcy, (23) 

which are sufficient to calculate all other quantities. The corresponding expression for Ri^, expanded at the leading 
power of Sz , is 



Ru 



2Cj;C5SyS^ C^C^S^ -\- si^5y 2s^Sy -\- SxSy 



+ 0(sz). (24) 
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In this pattern, to match the experimental results, the required values are {9x = 34'^, 9y = 42*^, 5 = 92.755*^, 6z = 
5"). For these inputs we obtain mi/ma = 1.050052527, ma/mg = 1.048504453, p = 87.72", a = 95° and 
Ri, = 0.033. The mass fitted from the observed Am^^j is = 0.156 eV. Then the derived values for the other 

remaining parameters are Am||_^ = 2.5 x 10"^ eV^, Mp = 0.164 eV, Mpp = 0.160 cV and S = 0.48 cV. There is 
no tuning required for the mixing angles 9x and 9y, to assure their consistency with the relation 0.025 < < 0.05 
as is shown in Fig. 2(a). In this pattern the numerically estimated mass matrix M is 



M 



ma 



Pattern D2: A/^ 

equation (20) are 



/ -1.0268 + 0.000489 i 0.027333 - 0.00049 i 0.21400 + 0.00037 i \ 

1.0268 - 0.000513 i -0.02737 + 0.000517i (25) 
-0.02737 + 0.000517 i 0.99942 - 0.000507 i j 

0. In this pattern the required quantities a's and Vs as given by 



0.027333- 0.00049 i 
V 0.21400 + 0.00037 i 



Mrr = 0, Mee + 



ai 
bi 



iS\2 



^x^z ~^ ( ^x^y^z 



-Z(5\2 
Sx 6 j , 



02 

b2 



S j.C^ I ^ S X ^y ^ z 



x^z I V ^x^y^z 
2^2 



Cx.Su^ ) 7 
iS\2 



^X O y 
Cx Cy G 



as 



2,22 



63 = sl 



2 2 



(26) 



which are sufficient to calculate all other quantities. 

Using Sz as a small parameter, we expand in terms of its powers and keep only leading terms. The analytical 
approximate formula for the mass ratio are: 



mi 
ma 



m2 
7713 



'"2x 
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The corresponding expressions for p and a are 



1 

- — tan 
2 



-1 



Oisz 



0{Sz 



while the corresponding expressions for Mpp and Mp are 

2Cy Sy^x 



1 

- tan 
2 



C5C,(l + 2c|c2 -24 -c2) 



2c^C5Si 



1 



0{sl 



0{Sz 



(27) 



(28) 



(29) 



and 



Mpp = ma + 0{sz), 



M, 



1 - Hshl 



0(s.). 



(30) 



In this pattern, to match the experimental results, the required values are {9x = 34°, 9y = 44.65°, (5 = 
90°, 9z = 5"). For these inputs we obtain mi/mg = 1.015274786, ma/mg = 1.014764820, p = 89.99°, a = 89.86° 
and R^ ~ 0.035. The mass mg fitted from the observed Am^^j is mg = 0.28 eV. Then the derived values for the 

other remaining parameters are Am^^^^ = 2.4 x 10"^ eV^, Mp = 0.28 eV, Mpp = 0.28 eV and E = 0.84 eV. The 
angle 9y is highly constrianed around 45° for a reasonable choice of 9^ as can be seen from Fig. 2 (b). In this 
pattern the numerically estimated mass matrix M is 



M = mg 



-0.99982 - 
0.12446 
0.12299 



0.0017117i 

- 0.00027 i 

- 0.00004 i 



0.12446- 0.00027 i 
0.99981 - 0.001714i 
-0.01515 + 0.0017276i 



0.12299 + 0.00004i 
-0.01515 + 0.0017276 i 
0.99982- 0.001714 i 



(31) 
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Figure 2: The available parameter space for the cases Di (a) with 5 = 92.7° and D2 with 5 = 90" (b) respectively. The dotted 
region represents the points {9x, 9y) at which Ri, is constrained by 0.025 < Rv < 0.05. The range from 29*^ to 57" is spanned for both 
angles 9x and 8y. 



Pattern D3: Af^^ 
equation (20) are 



Mrr = 0, Mfj^e + M^-t = 0. lu tliis pattern the required quantities a's and 6's as given by 



2 2 , 2 -2i(5 
L^rfO^ I l) ™ C- ^ t-iS 



2 -2i(5 



2 2 



as 



5\2 



5l ~ ( CxSySz SxCyC J^X^Z \~'^X^y'^ Z ~r ^a;OyC 

62 — ( ^x^y^z H~ C-xCyC )SxCz ~h ( SxCySz CxSyC ) , 
^3 — ^y^z^z ^y^z^ 



(32) 



which arc sufficient to calculate all other quantities. The corresponding expression for R^, expanded at the leading 
power of Sz, is 

2csCxCySx - Aslsl + 2sl + 2si - 1 



-2csCxCyS% + 2s^s' 



y"x 



+ 0(.s,). 



(33) 



In this pattern, to match the experimental results, the required values are {6x = 34", Oy = 42", 5 = 272.9°, Oz — 5"). 
For these inputs we obtain mi/mg = 1.008141341, ms/mg = 1.007871889, p = 177.55°, cr 5.32° and R^ = 0.034. 
The mass fitted from the observed Am^^^ is TO3 ~ 0.38 eV. Then the derived values for the other remaining 

parameters arc Am^^^ = 2.4 x 10"^ eV^, Mp ^ 0.38 eV, Mpp = 0.38 cV and S = 1.15 eV. In this pattern the 
numerically estimated mass matrix M is 



M = TO3 



/ 1.0001 - 0.000939 i -0.093575 - O.OOOOH 0.084241 + O.OOOIH \ 

-0.093575- O.OOOOH -0.09358 + 0.000006 i 0.99567 - O.OOOOOH 
\ 0.084241 + 0.00011 1 0.99567 - O.OOOOOH 0.09355 - 0.000011 i / 



(34) 



4 Normal hierarchy models 

Pattern Ni: Me^ + il/^r = 0,Mee + A/^t = 0. In this pattern the required quantities a's and 5's as given by 
equation (20) are 



ai 

0.2 

as 



= CxCz{- 
— ^xCz{~ 



-'X'^y'^z 



-'x^y^ 



) + (- 



~'X'-^y'=>z 



cCyG )( C^xC-ySz H~ ' 



CxCyG ) ( ^x^y^z 
bi ^x^z ( ^x^y^z SxCyG ' )( CxC-ySz 



SzSyCz 



JX'^y^z i 

■f SyC^Cy 



■^x ^y ^ 



^x '^y^ 



)(- 



=*x^y'^ z 



^X'-'y^ 



-iS\ 



7 



i5\ 



-iS\ 



■y'-> z ^x "^y^ 



63 = SyClCy, 



(35) 



which arc sufficient to caicuiate aii other quantities. 

Using Sz as a smaii parameter, we expand in terms of its powers and fceep only the leading terms. The 
analytical approximate formulae for the mass ratios are: 



mi 



I ^SxCxCSCy Cy 



44 , ^ rn2 

ma 



^Zc^CyCxSx ~t~ S^Sy 1 



+ 0(sz) (36) 



where 



-^ni — 4C(5S^Ca;Cy5y Q^SyCySxCxC§ 4:CsSyS^Cx ~\~ ^CgSySxCx H~ "^^y^y^S^x ^^V^V^x "^^V^y^S^x 

_L2r s - s^^s"* - Ss-^s^ + + s^'s^ + Ss^s^ - - 



(37) 



The corresponding expression for p and cr are 



1 

« — tan 
2 



+0(s^), (T«-tan — 



+ 0(s.), 



(38) 



where 



-^pl ~ ^x^si, ^^S^y^x^y ~l~ ^x^y^y ^~ ^^x^y^5^y ^~ '^a^ ^^S^y^x '^^S^y^x^x ~^ '^^S^y^x^x^y^ i 

Np2 = (4S5Q - 3Q)c2;S2/Cy4 + s4 + (1 - 2s^)(sy4 - s^s^ + 4*^) + C3;S^Q + csCyC^SxSy - Cysl, 

/22 2 32 322 \ 

Nal ^Sy^^x^y^S^y^X H" ^CxCySyC§ ^Cj.SyCyC§ ~ 'ZCgCxSy H~" ^C^C^Sy Cr^Sx H~" SyC^Sx.Cy CySySx)') 

Na2 = CgSx + Ca,CyS^ + (1 - 2S5)(S^C^S^ + C^Sj,s| - C^Sj^) - Qsj^ + (-4s^C5 + 3c5)Sy4s + - 2cs)SySxCy, 

(39) 



while the corresponding expressions for R^, Mpp and Mp are 



i?2 



+ 0(s,) 



(40) 



where 



where 



Ri = -2sxCxCyCsSy + 2s,yS^ - s,y 

323 24 4 '2. '2i "2 

R2 = 28-1 CxCyCSSy 4:SxCxCSSy + 2SxCxCSSy + 4SyCyS^S3, 2SyCyS^ ^^V^V^S^X ^SyCyS^. 

-344 + 4 + - si (41) 

M^^«m3^yX + o(,,), (42) 



and 



M- 



4 4 4_ o„4 2 _ 4 
°y°x ' '~'°x°y °x 



- s^ysl + [Aci - 2)cySysl + (2 - Aci)cySysl ^ + s . 

„3, 



AiC^SxC-xf^ySy 4" AiC^SyS j^Cx 2c§SySxC-x ~\~ 2SyCySxCxCS ^ 



Mj3 msCx] 



\sisl-2slsl , ^ 
' M2 



(43) 
(44) 
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where 



Mo 



2csSySxCx + 4 + (2 - ^c1)cySysl + 3s? 



-\-4:CSS.j.CxCySy 4:CSSyS,j,C'^ 



2 _ 2 _ 4 4 



2slCySxCxCS + (4C^ - 2)CySy 



Ss^Sy 



4 2 
^y^x 



(45) 



In this pattern, the required input values in order to match the experimental results are (O^ — 33'^,9y — 
42°, 6 138.1", 9^ ^ 5"). For these inputs we obtain mi/mg = 0.6563621852, 7712/^3 = 0.6416046770, p = 79.56°, 
a ~ 148.87° and i?^ — 0.033. The mass TO3 fitted from the observed Ato^^j^ is = 0.064 eV. Then the derived 
values for the other remaining parameters are ^'^i^^^jj^ = 2.4 x 10~^ eV^, M/3 = 0.042 eV, Mpp = 0.021 eV and 
E = 0.15 eV. As it is evident from Fig. 3 (a), the available parameter space permits the choice of 6x and 9y in 
the acceptable range without tuning. In this pattern the numerically estimated mass matrix M is 
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Figure 3: The available parameter space for the cases Ni (a) with 5 = 138.1'^and I4 (b) with S = 3l" respectively. The dotted 
region represents the set of points (^9x , Oy ) at which Ri, is constrained by 0.025 < < 0.05. The range from 29'' to 57° is spanned 
for both angles 6x and 9y. 



-0.33257- 0.00386 i -0.33256 - 0.003873 i 0.45632 + 0.003944 i 
M = m3 ( -0.33256- 0.003873 i 0.67411 - 0.00379H 0.33257+ 0.003870 i 
0.45632 + 0.003944 i 0.33257 + 0.003870 i 0.64681 - 0.003947 i 



(46) 



Pattern N2: M^^ + M^r — 0, M^e + ~ 0. In this pattern the required quantities a's and &'s as given by 
equation (20) are 



ai = 

as = 
bi = 
b2 = 
63- 



( ^x^y^z ^xCyG ' ^CxCz ~\~ ( CxSyS^ SxCyG )( CxCyS^ "4^ 



)( ^x^y^z 



( ^x^y^z ^~ 
3 — SySzCz ~\~ CyC^ 



-iS 



( ^x^y^z ^~ CxCyG ' '^SxCz ~f" ( '^x^y^z CxSyC ) 



^^x^z H~ ( ^x^y^z ~l~ ^x^y^ ' ) 5 



(47) 



which are sufficient to calculate all other quantities. The corresponding expression for R^, expanded at the leading 
power of s^, is 



^CxSj,CySyCS 2CxS j,CySyC§ ASyCxS^CS + QSyCx^xCS 2SyCxSxCS ^^y^y^x ^~ ^^y^y^x ^^y^x ~^ ^y^x 



D 



N2 



+ 0(5.). 

(48) 
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where 



+Qslcysl - 44c,4 - + - 44 + 444 - 6^4 + 44. (49) 

In this pattern, matching the experimental results, we find that the required input values are {Ox = 41. 5*^,^1, = 
47°, 5 = 195°, 9^ = 5°). For these inputs we obtain mi/mz = 0.7674926621, ma/mg = 0.7584041167, p = 96.62°, 
(7 = 8.4° and Ry ~ 0.033. The mass fitted from the observed Arrig^j is — 0.076 eV. Then the derived 

values for the other remaining parameters are ^'m-^i^ = 2.3 x 10"'^ eV^, Mp = 0.058 eV, Mpp = 0.007 eV and 
S = 0.19 eV. In this model 6x turns out to be a little bit out of the allowed range. In this pattern the numerically 
estimated mass matrix M is 

/ -0.092 - 0.0024 i -0.4444 + 0.000018 i 0.617 + 0.00029 i \ 
M = m3 -0.444 + 0.000018 i 0.639 + 0.000014 i 0.444 - 0.000018 i (50) 
\ 0.617 + 0.00029 i 0.444 - 0.000018 i 0.444 - 0.000018 i / 



5 Inverted hierarchy models 



Pattern Ii: Alee + ^-tfifi = 0, Alee + M^j^r = 0. The a's and 6's quantities can be inferred from the corresponding 
ones of cases Di and Ni. The corresponding expression for i?^, expanded at the leading power of Sz, is 



2CySy ACy 



q3 „4 



Dn 



+ 0{sz), 



(51) 



where 



Dn 



-l + (4 + 4c2)c,44 



8c,,s,4 + 44, + 24 + (-3 + 4c^)44 + (-4cJ + 
44, + ScySysl + (-8 + 4c^)44 + (-4c^ + 2)44 + (-2 - 4c^)c^44 



)44 



'i.CySy 



(52) 



In this pattern, to match the experimental results, the required input values are {Ox — 43°, 6y — 44°, 5 ~ 26°, Oz = 
5°). For these inputs we obtain mi/ms = 4.230294841, mz/mg = 4.298049784, p = 103.57°, a = 15.37° and 
Ru = 0.033. The mass TO3 fitted from the observed Am^^j is 777.3 = 0.012 eV. Then the derived values for the 

other remaining parameters are A?77|^j^ = 2.3 x 10~^ eV^, Mp = 0.05 eV, Mpp — 0.003 eV and S = 0.11 eV. In 
this pattern the numerically estimated mass matrix M is 



-0.2982 + 0.0109 i 3.1198 + 0.17560i -2.8549 - 0.17093 i 
M = ?773 ( 3.1198 + 0.17560 i 0.29818 - 0.01085 i 0.2982 - 0.01090 i 
-2.8549- 0.17093 i 0.2982 - 0.01090 i 1.0592 + 0.03131 i 



(53) 



Pattern I2: Mge + A/^^ = 0, M^e + -^^rr = 0. The a's and 6's quantities can be inferred from the corresponding 
ones of cases Di and N2. The corresponding expression for R^, expanded at the leading power of s^, is 



Rv 



-6s;, 



CyC^Sy 



34 



24 



0{Sz), 



(54) 



where 



D12 



-1 



(44 - 10)44 



6st 



+ (44 - 5)44 + (-44 + 12)4s 



4 

^CyCSSy 

2 _ 
x"y 



IcxCyCssl + (-44 + 2)sysl 



Ast + 2st, + 4s: + 2sxCxC. 



'yCSSy 



(55) 
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In this pattern, to match the experimental results, the required values are (6*^ = 34", 6y = 43°, S = 118°, 9z = 5"). 
For these inputs we obtain mi/ms = 2.062421918, mz/mg = 2.090403819, p = 77.05°, a = 132.41° and i?^ = 
0.034. The mass fitted from the observed Am^^j is ~ 0.026 cV. Then the derived values for the other 
remaining parameters are Am^^^ = 2.2 x 10~^ eV^, Mp = 0.054 eV, Mpp = 0.034 eV and E = 0.135 cV. In this 
pattern the numerically estimated mass matrix M is 



-1.3165- 0.03164i -1.0155 - 0.02210i 1.2241 + 0.02439 i 
M = m3( -1.0155- 0.02210 i 1.3165 + 0.03163 i -0.17371 - 0.02687i 
1.2241 + 0.02439 i -0.17371 - 0.02687i 1.0155 + 0.02213 i 



(56) 



Pattern I3: Me 



-M„, 



0, Muu + Mtt = 0. The a's and 6's quantities can be inferred from the corresponding 



ones of cases D3 and Ni . The corresponding expression for , expanded at the leading power of Sz , is 



SyCSSx( 



SyCSS^Cx 



(57) 



In this pattern, to match the experimental results, the required values are {d^ = 34°, 9y = 42°, S = 304°, 6z = 5°). 
For these inputs we obtain mi/mg = 1.676711363, mz/mg = 1.694189580, p = 166.78°, a = 47.34° and R^ = 
0.032. The mass 7713 fitted from the observed Am^^j is ~ 0.037 eV. Then the derived values for the other 
remaining parameters are Am^^^^ = 2.4 x 10"^ eV^, Mp = 0.062 eV, Mpp 0.036 eV and E = 0.160 eV. In this 
pattern the numerically estimated mass matrix M is 



0.98880 + 0.01487 i 
M = m3\ -1.0068- 0.00011 i 
0.90785 - 0.00166 i 



-1.0068 -O.OOOIH 
0.01346- 0.00008 i 
1.0068 + 0.00010 i 



0.90785- 0.00166 i 
1.0068 + 0.00010i 
-0.01346 + 0.00010 i 



(58) 



Pattern I4: M^^ 

equation (20) are 



A'Itt ~ 0, Mee + Mfj^r = 0. In this pattern the required quantities a's and 6's as given by 



ai 



i25 



0-3 



(59) 



where the coefficients 6's have the values as given in equation (35). 

Using Sz as a small parameter, expanding in terms of its power and keeping only leading terms, we have the 
analytical approximate formulae for the mass ratios: 



nil 

TO2 



1 2C^ Ac'^SyCy + AC^SyCy + 4:C'^Cy 4C^C^ 



-^x'^y^y^S 



+ 4+8cl4SyCy+0{Sz) 



~ C2iA/4c2 - AcisyCy - 44 + 4clsyCy " 8clcl + 8c2c4 + 44cg - 444 - SclsyCycj + 4 + ScjcisyCy + 0{Sz 



The corresponding expression for p and a are 



(60) 



1 

— tan 
2 



sls2S 



SIC2S - S2yCl 



+ 0(s.). 



1 

a K. — tan 
2 



-cls25 



S2ySl - CIC2S 



+ 0{Sz), 



(61) 



While the corresponding expression for R^, Mpp and Mp are 



1 + 44 



44 



^x^y ^Sj. 



(_8c2 + 4)s,c,4 - 2sl + 34 + (-4 + 8cj)syCysl + 4^4 _ 4,4,2 



+ Oisz), 



(62) 
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M/3/3 « ma S2y + 0{z), Mp « ms Cj^-VA/S, 



(63) 



where 



A/3 %C^^clSyCy - 



4cf, 



H (64) 



In this pattern, no need for tuning to match the experimental results as it is clear from the parameter space 
in Fig. 3(b). The required values are {dx — 34",6'j, = 42", (5 = 31°, 6*2 — 5°). For these inputs we obtain 
mi/ma = 1.601974769, ma/mg = 1.585960893, p 76.95", a = 134.67" and = 0.034. The mass TO3 fitted 
from the observed Am^^j^ is 7713 = 0.039 eV. Then the derived values for the other remaining parameters are 

^™atm = 2.4 x lO'^ eV^ Mp = 0.063 eV, Mpp = 0.039 eV and E = 0.165 eV. In this pattern the numerically 
estimated mass matrix M is 



-0.97926- 0.01129 i 0.12874 - 0.92428 i 0.11713 + 0.83357i 
Af = m3 ( 0.12874- 0.92428 i -0.10446 + 0.10832 i 0.97931 + 0.01130 i 
0.11713 + 0.83357 i 0.97931 + 0.01130 i 0.10442 - 0.10834 j 



(65) 



Pattern I5: M^e + M^t 0, A/^, 
of cases D2, while 6's are given as 



Mfj^T = 0. The a's quantities can be inferred from the corresponding ones 



hi = CxCz{CxCz - CxCySz + SxSyC * ), b2 ^ ~SxCz(-SxCz + SxCySz + CxSyC *''), bs, = Sz{Sz + CyCz). (66) 

The corresponding expression for R^, expanded at the leading power of Sz, is 



Ry 



^ 2Cj,c!j. + 2Sy ASyCx) 



Dl5 



OiSz), 



(67) 



where 



Di5 = - 



l + 4cl + 2cl + Aclclcj ~ -icjclclsy + 4cl4c1sy + ^c^Sy - e^c^s^ - %CySycl - Sclcl 

" 3c;)4 - 44 - ^.c^yclcl + ^0^x4 - 4cl4cj + 844 + 8cySyc2 + Ac^cl + Acl4sy 



'2CySy 



(68) 



In this pattern, there is no tuning for the angles 9x to match the experimental results. The required values are {dx = 
44.73", 6*^ = 42°, (5 = 90", 6*^ = 5°). For these inputs we obtain mi/ms = 4.137879853, m2/m3 = 4.076054427, 
p = 51.24", cr = 136.24" and i?,^ = 0.033. The mass nis fitted from the observed Am^^j is TO3 = 0.012 eV. 

Then the derived values for the other remaining parameters are Arri^^^ = 2.5 x 10^^ eV^, Mp = 0.051 eV, 
Mpp = 0.004 eV and S = 0.115 eV. In this pattern the numerically estimated mass matrix M is 



-0.35325 + 0.0220 i -2.9235 - 0.39769 i 2.7916 + 0.35549 i 
Af = m3 ( -2.9235- 0.39769 i 0.98990 + 0.0765 i 0.35326 - 0.0220 i 
2.7916 + 0.35549 i 0.35326 - 0.0220 i 0.35325 - 0.02201 i 



(69) 



Pattern le: A/ee+A/^^ = 0, A/^^+A/t-t- = 0. The a's and 6's quantities can be inferred from the corresponding 
ones of cases D2 and D3. The corresponding expression for R^, expanded at the leading power of Sz, is 



Ru 



244- 



st. + : 



x"y 



''y^x 



„2„2„4 
^S^y'^x 



„4 „2 4_ 2 
^x^S + ^x 



■^x ^x'^y 



OiSz). 



(70) 



In this pattern, no tuning is needed to match the experimental results. The required values are {9x — 42",0.y = 
42", 5 = 170", 6*;, = 5"). For these inputs we obtain mi/mg = 8.699398197, 1112/1113 = 8.840977186, p = 174.85", 
cr = 83.86" and i?^ = 0.032. The mass 7713 fitted from the observed Ato^^j is TO3 = 0.006 eV. Then the derived 
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values for the other remaining parameters are ^"t,|^|.^ = 2.4 x 10 eV^, Mp 
S ~ 0.105 cV. In this pattern the numericaUy estimated mass matrix M is 



0.049 eV, Mpp = 0.005 eV and 



0.8601 - 0.01705 i 6.4626 -0.1446H -5.8026 + 0.13222 i 
M = m3( 6.4626- 0.14461 i -0.8601 + 0.01705 i 0.9140 + 0.00166 i 
-5.8026 + 0.13222 i 0.9140 + 0.00166 i 0.8602 - 0.01705 i 



(71) 



Pattern I7: Mee + Mrr = 0, Mefi + Mrr = 0. The a's and b's quantities can be inferred from the corresponding 
ones of cases D2 and D3. The corresponding expression for R^, expanded at the leading power of Sz, is 



Cy{2csSxCxCy - 1 + 2cl) 



(72) 



where 



D 



77 



^222 



-'y 



„2 2 



I2csclcls^ - 2c2 + ^csclcyS^ + ic^c. 



4 



^44 



,3 

y 



y ^ ^44"^^ 

^^x^y^5^x 



(73) 



In this pattern, no tuning for the angles Ox to match the experimental results. The required values are {Ox ~ 
31", Oy = AA°,S = 147°, 6*^ = 5°). For these inputs we obtain mi/mg = 1.475002097, mz/mg = 1.462135796, 



0.046 eV. 



p = 81.84°, a — 154.74° and R,y = 0.033. The mass ma fitted from the observed Am^^j is TO3 
Then the derived values for the other remaining parameters are Arri^^^ = 2.5 x 10^^ eV^, = 0.067 eV, 
Mf^l^ = 0.036 eV and S = 0.180 eV. In this pattern the numerically estimated mass matrix M is 



-0.77989 + 0.00498 i -0.77991 + 0.00497 i 0.96962 - 0.00542 i 
M^ni-i ( -0.77991 + 0.00497 i 0.98867 - 0.00664 i 0.10578 + 0.00582 i 
0.96962- 0.00542 i 0.10578 + 0.00582 i 0.77989 - 0.00498 i 



(74) 



Pattern !§: M/^if^+Mrr = 0, Mee+^J^rr = 0. The a's and &'s parameters can be inferred from the corresponding 
ones of cases D2 and D3. The corresponding expression for R^, expanded at the leading power of s^, is 



Ri> 



4sf 



Asi 



- 8^4 



4'SySx 



8^2452 + ; 



^iS^gS^Sx 



■As 



2 „2 
x^y 



8s^Sy 



2 + 4s? 



''x^'S 



Asi 



x^S 



0(s. 



(75) 



In this pattern, no need for tuning for the angles Ox to match the experimental results. The acceptable values 
arc {Ox = 42.9°, 6lj^ = 42°, (5 = 33°, = 5°). For these inputs we obtain mi/ms = 10.88384938, ma/mg = 
11.06236809, p = 18.5", a = 110.84° and R^ = 0.032. The mass mg fitted from the observed Am^^j is mg = 

0.004 eV. Then the derived values for the other remaining parameters arc ^™^tm ~ ^'^ ^ 10"'^ eV^, A/^ = 
0.049 eV, Mg/3 = 0.004 eV and S = 0.103 eV. In this pattern the numerically estimated mass matrix M is 



M 



mg 



0.8378- 
-8.0358 
7.2547 4 



- 0.1050 i 

- 0.90178 i 
0.79963 i 



-8.0358 
0.8375- 
1.0923- 



- 0.90178 i 

- 0.1049 i 
-0.0118i 



7.2547 + 0.79963 i 
1.0923 + 0.0118 i 
-0.8376- 0.1048 i 



(76) 



The remaining two models of Mn + Afgg = 0, M12 + M23 = and Afn + M23 = 0, M21 + Afgg = fail to be 
consistent with the experimental data for any reasonable choice of the two mixing angle Ox and Oy . 
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Model 


Trace conditions 


Status 




Oy 


Oz 


r 





7711 


?Tl2 


P 


(T 


ma 


Mp 




2j 


A 1 


Di 


(11, 22), (12, 23) 


allowed 


34 


42 


5 


92.755 


0.033 


1.050 


1.049 


87.72 


95 


0.156 


0.164 


0.160 


0.484 


0.0025 


D2 


(11, 33), (11, 22) 


allowed 


34 


44.65 


5 


90 


0.035 


1.0152 


1.0147 


89.99 


89.86 


0.276 


0.281 


0.277 


0.84 


0.0024 


D3 


(22, 33), (21, 33) 


allowed 


34 


42 


5 


272.9 


0.034 


1.0081 


1.0079 


177.55 


5.32 


0.38 


0.38 


0.38 


1.15 


0.0024 


Ni 


(12, 23), (11, 23) 


allowed 


33 


42 


5 


138.1 


0.033 


0.656 


0.642 


79.56 


148.87 


0.064 


0.042 


0.021 


0.148 


0.0024 


N2 


(12, 23), (12, 33) 


disallowed 


41.5 


47 


5 


195 


0.033 


0.767 


0.758 


96.62 


8.4 


0.076 


0.058 


0.007 


0.191 


0.0023 


Ii 


(11, 22), (11, 23) 


disallowed 


43 


44 


5 


26 


0.033 


4.230 


4.298 


103.57 


15.37 


0.012 


0.05 


0.003 


0.111 


0.0023 


I2 


(11, 22), (21, 33) 


allowed 


34 


43 


5 


118 


0.034 


2.0602 


2.0904 


77.05 


132.41 


0.026 


0.054 


0.034 


0.135 


0.0022 


I3 


(12, 23), (22, 33) 


allowed 


34 


42 


5 


304 


0.032 


1.677 


1.694 


166.78 


47.34 


0.037 


0.062 


0.036 


0.160 


0.0024 


I4 


(22, 33), (11, 23) 


allowed 


34 


42 


5 


31 


0.034 


1.6019 


1.5859 


76.95 


134.67 


0.039 


0.063 


0.039 


0.165 


0.0024 


I5 


(11, 33), (23, 11) 


disallowed 


44.73 


42 


5 


90 


0.033 


4.137 


4.076 


51.24 


136.24 


0.012 


0.051 


0.004 


0.115 


0.0025 


le 


(11, 22), (22, 33) 


disallowed 


42 


42 


5 


170 


0.032 


8.699 


8.841 


174.85 


83.86 


0.006 


0.050 


0.005 


0.105 


0.0024 


I7 


(11, 33), (21, 33) 


allowed 


31 


44 


5 


147 


0.033 


1.475 


1.462 


81.84 


154.74 


0.046 


0.067 


0.036 


0.180 


0.0025 


Is 


(11, 33), (22, 33) 


disallowed 


42.9 


42 


5 


33 


0.032 


10.88 


11.062 


18.5 


110.84 


0.004 


0.049 


0.004 


0.103 


0.0024 



Table 1: The 13 patterns for the two- vanishing traces. The trace corresponding to the index {ab, ij) is Mab + Mij = 0. All the angles are measured in degrees, masses in 
eV and AmL„ in eV^ 



6 Conclusion and discussion 



In this work several patterns of Majorana neutrino mass matrix, consistent with the present available observed 
data, are derived. The models are based on textures possessing two 2x2 sub-matrices with vanishing trace. The 
new proposed texture can be considered as a non trivial generalization of the zero-texture as explained in the 
introduction. 

In our work we have thirteen possible acceptable patterns for Majorona mass matrix out of fifteen ones. 
The resulting models fall into three distinct classes namely, degenerate case (Di • • - Da), normal hierarchy case 
(Ni andN2), and inverted hierarchy case (Ii • • - Is). The numerical results of our study is summarized in Ta- 
ble (l)and Table (2) for a quick reference. 

Our numerical study reveals that there are eight models (Di, D2, D3, Ni, I2, 13, 14 andlr), for which the 
mixing angles {9x, Oy andOz) can be adjusted to fall into the acceptable range given by equation (13). In the 
remaining five models, (N2, Ii, I5, le and Is), Ox falls in the range 41.5° < 0^ < 45°, which is out of the acceptable 
range. These kinds of models can be considered as empirically ruled out. An avenue for curing these models could 
be provided by a small perturbation over the adopted textures. In our subsequent discussion we only focus on 
the successful models namely (Di, D2, D3, Ni, I2, 13, 14 andir). 

The numerical study points out that Dirac-type phase tends to be around ^ or ^ in the degenerate 
case, while in other cases no general trend could be observed. Concerning a possible relation between Dirac and 
Majorana phases that could be revealed by numerical study, we find the two phases p and d almost satisfying the 
relation p « | in the normal hierarchy case, while for the degenerate case there is the relation p ^ S or S — ^, 
and for the inverted hierarchy case p « | is obeyed except for I4 where p « 2 (5 is satisfied. 

Another possible relation between a and S could be easily recognized. In the normal hierarchy and degenerate 
case, the relation a Ri S or S ~ ^ is satisfied. The inverted hierarchy cases have no specific general relation which 
is obeyed. 

The non oscillation parameters M^,M^^andS are consistent with the bounds given in equation (16). In all 
the successful models, M3, Mpp and ms have the same order of magnitude. The mass sum parameter is always 
constrained to be S < 1.15 eV which is safe with the cosmological bound in equation (16). 

All successful models are found to be still consistent with experimental data in the limit of vanishing 6z while 
keeping 0^ ~ 34° and 6y w 42° constants. The same thing still holds, when 9^ is stretched to its upper bound 
(10°). In these limits, little changes take place for the other parameters. 

Regarding the hierarchical structure of the mass matrices, as it is evident from table. 2, all successful models 
have clear hierarchical structure except the mass matrix of model Ni whose elements have all the same order 
of magnitude. These hierarchical properties are restricted to the real parts, but for the imaginary part they are 
always very small in comparison with the real ones with the exception of model I4. 

Final remark, related to when we restrict the study to the parameter space {O^ ~ 34'-', 6y « 42° and 6^ « 5°), 
while varying 6 under the condition 0.025 < Ru < 0.049, we find that all successful models turn out to be tightly 
constrained in order to have a quasi degenerate spectrum, (mi ~ TO2), for mi and m2 and no strong hierarchy 
between mi ^ m2 and ma can occur. This can be considered as a general prediction for these class of models. 
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Model 


Trace conditions 


M 












-1.0268 + 0.000489 i 


0.027333 - 


0.00049 i 


0.21400 + 0.00037 i 






Di 


(11, 22), (12, 23) 




0.027333- 0.00049 i 


1.0268 - 0.000513 i 


-0.02737 + 0.000517 


i 














0.21400 + 0.00037 i - 


-0.02737 + 


0.000517 i 


0.99942- 0.000507 i 














-0.99982 + 0.0017117 i 


0.12446- 


- 0.00027 i 


0.12299 + 0.00004 


i 




Da 


(11, 33), (11, 22) 






0.12446- 0.00027 i 
0.12299 + 0.00004 i 


0.99981 - 
-0.01515 + 


0.001714 i 
0.0017276 i 


-0.01515 + 0.0017276 
0.99982- 0.001714 i 


i 














1 nnm — nnnq'^q 7 


-0.093575- 0.00001 i 


084941 4- 0001 1 






Da 


(22, 33), (21, 33) 






-0.093575 - O.OOOOH 


u.uyooo 


1 n ooooofi i 


0.99567- 0.000001 


i 














0.084241 + 0.00011 7 


0.99567 - 


C\ A A A A A 1 ■ 

- U. 000001 % 


0.09355- 0.000011 


i 














— n 339^7 — 00386 7 




U.UUOO ( i 


45639 -1- 003944 






Ni 


(12, 23), (11, 23) 






-0.33256- 0.003873 7 


n fi741 1 - 

U .U ( 4:± ± 


- OO^^TQl ?■ 


0.33257 + 0.003870 


i 














0.45632 + 0.003944 7 




- U.OOooyO I 


0.64681 - 0.003947 


i 














—0 092 — 0024? 


— n 4444 - 

U. 4:4:4:4: 


4- 00001 S 7 


61 7 00029 7 






Na 


(12,23), (12,33) 




7773 


-0.444 + 0.000018 7 




00001 4 ? 


0.444- 0.000018 7 
















0.617 + 0.00029 7 


0.444 — 


A A A A A 1 " 

0.000018 2 


0.444-0.000018 7 
















n 9QS9 ^ n m no 1 












Ii 


(11, 22), (11, 23) 




777,3 


3.1198 + 0.175607 




- 01 OR^ 7 


0.2982 - 0.01090 7 
















-2.8549- 0.17093 7 


0.2982 — 


A A 1 A A A 

0.01090 2 


1.0592 + 0.031317 




















0991 7 








I2 


(11, 22), (21, 33) 


7713 




-1.0155- 0.022107 




0*^1 Pi'X 7 


-0.17371 - 0.02687 i 
















1.2241 + 0.02439 7 


—0.17371 - 


A A 0/^ ^ ■ 

- 0.02687 z 


1.0155 + 0.022137 
















u. jooou u.u±4:0 ( fc 


1 .uuuo 


0001 1 7 








I3 


(12,23), (22,33) 




7773 


-1.0068 - 0.00011 7 




OOOOJ^ 7 


1.0068 + 0.00010 7 
















0.90785 - 0.00166 7 


1.0058 + 


A A A A 1 A ' 

O.OOOIO7 


-0.01346 + 0.00010 7 


















n 1 9874 
u. izo ( 4 


n Q949S i 








I4 


(22, 33), (11, 23) 


7713 


0.12874-0.92428 7 


— 10446 


4- 1 0839 7 


0.97931 + 0.01130 i 
















0.11713 + 0.833577 


u.y /yoi 


+ U.UiioU 7 


0.10442 - 0.10834 7 
















—0 35325 + 0220 ? 


-2.9235- 


- 0.39769 i 


9 701 R 1 n '^5549 7 






I5 


(11, 33), (23, 11) 




77I3 


-2.9235- 0.39769 i 


0.98990 


+ 0.0765 i 


0.35326- 0.0220 i 
















2.7916 + 0.35549 i 


0.35326 


- 0.0220 i 


0.35325- 0.02201 i 
















0.8601 - 0.01705 i 


6.4626- 


0.14461 i 


-5.8026 + 0.13222 i 






le 


(11,22), (22,33) 




77I3 


6.4626- 0.14461 i 


-0.8601- 


h 0.01705 i 


0.9140 + 0.00166 i 
















-5.8026 + 0.13222 i 


0.9140 + 


0.00166 i 


0.8602- 0.01705 i 
















-0.77989 + 0.00498 i 


-0.77991 


+ 0.00497 i 


0.96962- 0.00542 i 






I7 


(11, 33), (21, 33) 


7713 


-0.77991 + 0.00497 i 


0.98867 


- 0.00664 i 


0.10578 + 0.00582 i 
















0.96962- 0.00542 i 


0.10578 


+ 0.00582 i 


0.77989- 0.00498 i 
















0.8378 + 0.1050 i 


-8.0358 


- 0.90178 i 


7.2547 + 0.79963 i 






Is 


(11, 33), (22, 33) 




7773 


-8.0358- 0.90178 i 


0.8375- 


f 0.1049 i 


1.0923 + 0.0118i 
















7.2547 + 0.79963 i 


1.0923- 


f 0.0118i 


-0.8376- 0.1048 7 







Table 2: The mass estimates for the 13 acceptable patterns. The trace corresponding to the index {ab,ij) 
is Mab + Mij = 0. 
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